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Abstract — Peristaltic transport of a conducting 
Jeffrey fluid in an inclined channel with heat transfer 
is studied under long wave length and low Reynolds 
number approximations. An exact solution is 
presented for the temperature field. The expressions 
for axial velocity, axial pressure gradient, Stream 
line and coefficient of heat transfer have been 
obtained analytically. The effect of various 
parameters on the flow characteristics are 
discussed with the help of graph. 
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I. INTRODUCTION 



Peristaltic pumping is a form of fluid transport which 
occurs in a biological system. The study of peristaltic 
motion has gained considerable interest because of its 
extensive role in transporting many physiological fluid in 
the body in various situations such as urine transport 
from the kidney to bladder, transport of the spermatozoa 
in the ducts afferents of male reproductive tract, the 
moment of chyme in the gastrointestinal tract, 
swallowing of food through esophagus and vasomotion 
of small blood vessels. Many mechanical devices have 
been designed on the peristaltic pumping to transport 
fluids without internal moving parts, for example, the 
transport of slurries, sensitive or corrosive fluids, 
sanitary fluid and noxious fluids in the nuclear industries. 
The literature on this topic is quite extensive. Mention 
may be made to some recent theoretical and 
experimental contributions Shapiro and Jaffrin [1], Hayat 
and Sajid [2], El Misery et al.[3],Yin and Fung[4,5], 
Pzrikidis [6], Shukia and gupta[7], Vjravelu et al.[8,9], 
Raju and Devanathan[10] in the field for Newtonian 
fluids. Such approximation is true in the ureter but it fails 
to give an adequate understanding of food mixing and 
chyme moment in the intestine, fuel slurries, flow of 
plasma, flow of mercury amalgams, and lubrication with 
heavy oils and greases etc. Also, the assumptions that 
most of the physiological fluid behave like Newtonian 
fluid is not true in reality. With all these facts in mind, it is 
clear that non Newtonian fluid plays an indispensable 
role in peristaltic flow problem. 



A number of researchers have been discussed 
the effect of magnetic field on peristaltic flow Mekheimer 
[11], Harikrishana and Subbareddyy [12], 
Suryanarayana Reddy et al. [13], Subba Reddy and 
Gangadhar [14], Elshahed and Haroun [15], 
Kothadapani [16, 17], Affi and Gad[118], Krishana 
kumari and Raman Murthy [19], Nadeem and Akbar 
[21], due to its applications in bio engineering and 
medical devices. Specially, magnetic wound or cancer 
treatment causes hyperthermia, bleeding reduction 
during surgeries, and targeted transport of drugs using 
magnetic particles as drug carries are few examples. 

The interaction of peristalsis in connection with 
heat transfer has also received some attention Nadeem 
and Akbar [21], Srinivas and Gayatri [22], Hayat and 
Heena [23], Hayat and AN [24], Hayat et al.[25] as it 
might be relevant in processes like hemodialysis and 
oxygenation. Rathod and Mahadev [26] have studied the 
effect of thickness of porous material on the peristaltic 
pumping of Jeffrey fluid with non-erodible porous ling 
wall. Effects of magnetic field and an endoscope on 
peristaltic motion is studied by Rathod and Asha [27]. 
Rathod and Pallavi [28] have discussed the influence of 
wall properties on MHD peristaltic transport of dusty 
fluid. 

The present article studies the effect of heat 
transfer on peristaltic MHD flow of Jeffrey fluid in an 
inclined channel. The governing equations of Jeffrey 
fluid in Cartesian co-ordinate have been modeled. The 
equations are simplified using long wave length and low 
Reynolds number approximations. The closed form of 
solution for velocity field and temperature are obtained. 
The influence of various parameters on the flow 
characteristics, the temperature and heat transfer 
coefficient are discussed through graph. 

2. FORMULATION OF THE PROBLEM 

Consider the peristaltic pumping of a Jeffrey 
fluid with heat transfer in an inclined channel of half- 
width 'a'. A longitudinal train of progressive sinusoidal 
wave takes place on the upper and lower wall of the 
channel. We further assume that the fluid is electrically 

conducting. A uniform magnetic field is applied in 

the transverse direction to the flow. The Reynolds 
number is taken small so that the induced magnetic field 
is neglected. For simplicity, we restrict our discussion to 
the half width of the channel as shown in figure .(1) 
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Figure. 1. Physical model 

The wall deformation is given by 

In 

H{x,t)=a + bCos — {X-ct) (1) 

A. 

where b is amplitude of the waves and A is the wave length. 
The constitutive equations for an incompressible Jeffrey fluid are 

T = -PI+S (2) 
S=^{7 + ^2Y) (3) 



where T and S are Cauchy stress tensor and extra stress tensor, p is the pressure, / is the identity tensor, // is 

the dynamic viscosity, \ is the ratio of relaxation to retardation times, is the retardation time, y is the shear 

rate and dots over the quantities denote differentiation. 

Under the assumptions that the channel length is an integral multiple of the wave length X and the pressure 
difference across the ends of the channel is a constant, the flow is inherently unsteady in the laboratory frame 

(X,7)and become steady in the wave frame which is moving with velocity 'c' along the wave. The 

transformation between these two frames is given by 

x = X-ct,y = Y,u = U,v = V (4) 

where U and V are velocity components in the laboratory frame and u and v are the velocity components in the 
wave frame. In the many physiological situations it is proved experimentally that the Reynolds number of the flow is 
very small. So, we assume that the wavelength is infinite, the flow is of Poiseuille type at each local cross-section. 
Introducing the non-dimensional quantities 

Ittx y u V ^ Ittu Ina'p 

A a C CO A JUCA 

t = ^,S = ^S,<l>=^-,G = P^,v = ^,R = '-^,Pr = ^, 
X idc a jjLc p V k 

where R is the Reynolds number, G is the gravitational parameter, £'c is the Eckert number and Pr is the Prandtl 
number. 

The equations governing the flow become 

du ^du 

— + S — = 0, (5) 

dx dy 
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and 



(6) 



(7) 



(8) 



1 



I 5y 1 + dy- 



The non dimensional boundary conditions are 

du . dO ^ . 
— = 0; — = 0 at y = Q 

dy dy 

u = —l ,0 = 1 at y = h 

Using long wavelength approximation and dropping terms of order S and higher. 
Eqs. (6) - (8) reduces to 



d 



1 du 



dy\\^X dy 

^ = 0 

dy 

1 d'e 



J 



- (w + 1) + GSina = 



dp 
dx 



Pr dy' 



+ Ec 



= 0 



where M = — — 



3. SOLUTION OF THE PROBLEM 



(9) 
(10) 



(11) 



(12) 



(13) 



Solving the Eq. (1 1 ) with the boundary conditions (9), we get 

P — Gsina 



u = 



M^coshiMh-Jl + I) 



cosh(iM^T+I)y) + 



^ -P + Gsina -M^^ 



V 



M 



J 



where P = 



dp 
dx 



(14) 



Integrating the equation (14) and using the conditions i// = 0 aty = 0 , we get the stream function as 

^ (-P + G sin a) sinh hMyll + A^ 



y/ =\{{y cosh[/;M7l + /I ]) (M' + P - G sin «) + 



V 



(15) 



(cosh W-s/1 + ;i + sinh /? JM" (1 + /I ) - 1 + Tank J 0- + ^)]} 



Substituting Eq. (14) into Eq. (13) subject to the boundary conditions (10), the temperature is 
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e = --^ {2M' (2M' - Br(p - G sin af (h^ - )(1 + X)) + (4M' + Br{p -Gsinaf) ^^^^ 



+ cosh(2/iM>/m) - Br{p - G sin af cosh(2Mj;Vm) sech(M/^^/fTX)' } 

where 5r = £'cPr is the Brinkman number. 
The coefficient of heat transfer is given by 

z = he . 

X y 

The volume flux 'q' through each cross section in the wave frame is given by 

h 



(17) 



q = judy (18) 

0 



1 {P-Gsina)Tanh{M^\ + X)h 



The expression for pressure gradient from Eq. (19) is given by 

dp _ M^/l + l((^)M' + (M' - Gsina )h) + Gsina.Tanh{M^\ + Xh) 



i^P - Gsina + M^^h (19) 



(20) 



dx Tanh{M^\ + Ah) - Mh^JT+I 

The instantaneous volume flow rate Q(x,t) in the laboratory frame between the central line and the wall is 

h 

Q(x, t) = j(u-\- \)dy = q + h (21) 

0 

Averaging the Eq. (21) over one period yields the time mean flow rate (time averaged flow rate) Q as 

1 r 

Q = -\^Qdt = q + \ (22) 
The pressure difference and frictional force across one wave length in an inclined channel is given by 



Ap={^^dx (23) 
^ Jo 



F=[^h-—dx (24) 
Jo 



'0 dx 



4. Results and Discussion 

To Study the behavior of the distribution of the 
axial velocity u, numerical calculations for several values 

of Jeffrey fluid parameter /I , Hartman number M , angle 
of inclination a , gravitational parameter G are carried 

out. Fig. 2(a) shows that an increase in X results in 
increase of velocity distribution. The effect of a on the 
velocity distribution can be seen through Fig. 2(b). It 
reveals that the axial velocity decreases with increasing 
a . Fig. 2(c) displays the effect of G on the velocity 
distribution for fixed values of other parameters. It is 
observed that the velocity increases with increase of G. 
The effect of the M on the velocity distribution is 
illustrated in Fig. 2(d). It is evident that, increase the 
value of M has a tendency to slow down the fluid motion 
and fluid moves like a block, which shows some sort of 
rigidity. This is because of the presence of the 
transverse magnetic field creates a resistive force 
similar to the drag force that acts in the opposite 
direction of the fluid motion, thus causing the velocity of 
the fluid to decrease. 



The effect of heat transfer on peristalsis is 
illustrated in Fig. 3. Fig. 3(a) is made to see the variation 

temperature 9 for various values of Brinkman number 

Br . It is observed that the temperature profiles are 

almost parabolic and the increase of 5r , the 
temperature distribution increases. From Fig. 3(b), it can 
be noticed that the temperature decreases with an 
increase of M. Figs. 3(c) and 3(d) are plotted to see the 
influence of Jeffrey fluid parameter and angle of 
inclination on the temperature distribution. It is observe 

that an increase in Jeffrey fluid parameter X affects the 
temperature profile in an opposite way to that of angle of 
inclination. 

Fig. 4 is plotted to study the effect of Jeffrey fluid 
parameter X , Hartman number M , angle of inclination 
a , gravitational parameter G on the pressure gradient 
Figs. 4(a) and 4(c) indicate that the pressure gradient 
increases with increase of M and G. From Figs. 4(b) and 
4(d), it is observe that pressure gradient decreases with 

increase of a and X . 
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Fig. 5 shows the behavior of heat transfer co- 
efficient Z. From this figure it is observed that due to 
peristalsis the heat transfer coefficient is in oscillatory 
behavior. The absolute value of heat transfer coefficient 

decreases with increasing aandM while it increases 
with increase of Br and X . 



5. TRAPPING 

Trapping is an interesting phenomenon in 
peristaltic motion. It is basically the formation of an 
internally circulating bolus of the fluid by closed stream 
lines. The effect of M on trapping can be seen in Fig. 6. 
We observed that the size of the bolus reduces with an 

increase in M. The effects of A and G on the stream 
lines are plotted in Figs. 7-8. It is observed that the size 
of trapping bolus increases with increasing A and G. 
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Figure. 2. The velocity distribution for = 0.4 , x = 0.25 ); other parameters are {a)a = —,G = 0.5, M = 2 , (b) 

M = 2,G = 0.5,^ = 1, (c)a = -,M = 2,A = l,(d) a = -,G = 0.5,A = l. 
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Figure. 4. Pressure gradient versus x for (a) ^ = 0.4 ,/l = 2,a = —,G = 3 . 
(b) G = 2,A = 2,M = 1,(^ = 0.4. (c)a = -, A = 2,M = !,(/) = 0.4. (6) a = -,G = 2,M = l(/) = 0.4 . 
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Figure. 5. Coefficient of heat transfer for (a) ^ = 0.6, M = 5, G = 0.2, A = 2, Br = 3, 

(b) a = -, G = 02,Pi = 2, Br = 3, (/) = 0.6. (c) a = -, G = 0.2, A = 2, M = 5, (/) = 0.6 

(d) a = -, G = 0.2, Br = 3,M = 5,(/) = 0.6. 
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Figure. 6. Streamlines for (a) M = 0.1.(b) M = 0.2.(c)M = 03. (d)M = 0.4 ; other 

parameters are <p = 0.4,0 = 30^ ,A = 4,G = 0.S . 
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c d 

Figure. 7. Streamlines for (a) /I = 2 (b) A = 3 {c)A = 4 {6) A = 5; other parameters are 

^ = 0.4, 0 = 3O\M = 0.4, G = 0.8 
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Figure. 8. Streamlines for different values of (a) G = 4 (b) G = 5 {c)G = 6 
(d) G = 7 ; Other parameters are = 0.4, 0 = 30\M = 0.4, A = 2. 
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